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A novel MOND effect in isolated high-acceleration systems 

Mordehai Milgrom 
DPP A, Weizmann Institute 

ABSTRACT 



I discuss a novel MOND effect that entails a correction to the dynamics of isolated 
mass systems even when they are deep in the Newtonian regime: systems whose extent 
R <C r M , where r M = (GM t / a ) 1 / 2 is the MOND radius, and M t the total mass. 
Interestingly, even if the MOND equations approach Newtonian dynamics arbitrarily 
fast at high accelerations, this correction decreases only as a power of R/r M . The effect 
, appears in formulations of MOND as modified gravity, governed by generalizations 

of the Poisson equation. The MOND correction to the potential is a quadrupole field 
Q_i! <p a ~ GQijr l r 3 , where r is the radius from the center of mass. In QUMOND (quasilinear 

MOND), Qij = —aQijr^j 5 , where Qij is the quadrupole moment of the system, and 
a > is a numerical factor that depends on the interpolating function. For example, 
the correction to the Newtonian force between two masses, m and M, a distance £ 
apart {£ <C r M ) is F a = 2a(£/r M ) 3 (mM) 2 M^ 3 a (attractive). Its strength relative to 
the Newtonian force is 2a(mM/M^)(a /g N ) 5/2 (g N = GM t /£ 2 ). For generic MOND 
f — \ theories, which approach Newtonian dynamics quickly for accelerations beyond o , the 

m 



predicted strength of the effect in the Solar system is rather much below present testing 
capabilities. In MOND theories that become Newtonian only beyond kcl , the effect is 
enhanced by k 2 . 

Subject headings: 



5_i , 1. introduction 



One would suppose (as had I for many years) that MOND dynamics always approach New- 
tonian dynamics, in regions of high accelerations, as precipitously as its field equations approach 
those of Newtonian dynamics. (For a recent review of MOND see Famaey & McGaugh 2012.) One 
would surmise that if the interpolating function that characterizes the theory, fi(x), approaches 
unity, very quickly, at high x, then the field, the forces, the dynamics in general, approach their 
Newtonian values as quickly, when the system accelerations become high compared with the MOND 
constant a . This may be so in some formulations (e.g., in 'modified-inertia' ones) and it is so for 
spherical systems, in general. But, I show here that this is not the general rule in formulations of 
MOND as an extension of the Poisson equation: the nonlinear Poisson formulation of Bekenstein 
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& Milgrom (1984), and the Quasi-linear 1 formulation (QUMOND) of Milgrom (2010a). No matter 
how vanishingly small 1 — \i is for high a; -namely how close the field equations themselves approach 
the Poisson equation at high accelerations-in aspherical systems, g/g N — 1 remains finite, and is of 
order (a /<?) 5//2 ; here g and g N are the MOND and Newtonian accelerations, respectively. 

Previously noted MOND corrections in systems with g S> a , are of several distinct types: 

(1) Corrections stemming from remaining departure of fi(x) from unity even for high acceler- 
ations (x 3> 1). Such corrections have been discussed extensively since the early days of MOND 
(e.g., Milgrom 1983, Sereno & Jetzer 2006, Milgrom 2009a). I do not discuss these here. To isolate 
away such effects, I shall assume here that fi(x) — 1 — > fast enough as x — > oo (how fast will 
become clear below). We can then put fi = 1 for x S> 1, to the desired accuracy. 

(2) Effects stemming from a MONDian background field in which the system is embedded, 
such as the galactic field in the context of a stellar system. This effect and its consequences for the 
Solar system have been considered in detail in Milgrom (2009a) and in Blanchet & Novak (2011), 
and these too are not considered in this paper, where I assume that the system is isolated. 

(3) Even in systems that are of high acceleration almost everywhere, there are, generically, 
small regions where the accelerations are smaller than a . These include, e.g., the very centers of 
stars and the regions around the points of zero gravity in many-body systems. Such near-zero- 
gravity regions in the Solar system have been discussed as possible sites for testing MOND using 
test-particles probes (Bekenstein & Magueijo 2006, Magueijo & Mozaffari 2012, Galianni & al. 
2011). Here I am not interested in the direct probing of such regions, but, rather, in the effect their 
presence has on the dynamics of the mass sources themselves. I shall show that the effects of these 
zero-gravity points are, generically, much smaller than the new effect I discuss here. 

In section 2, I describe the new effect qualitatively. In section 3, I calculate it in the framework 
of QUMOND. In section 4, I comment on the effect in the nonlinear-Poisson formulation. Section 
5 treats the zero-gravity points in the system. In section 6, I mainly discuss the effect in theories 
that restore Newtonian dynamics only at accelerations much larger than a . 

2. Qualitative explanation 

Consider an isolated distribution of nonrelativistic masses, p(r), of total mass M$, contained 
within a region of radius R, much smaller than the MOND radius of the system, r M = {M t G / (to) 1 / 2 . 
Also assume that the system is quasi-static, in the sense that it varies only on time scales much 
longer than r M /c. Apart for small regions around zero-gravity point-which I consider separately 
in section 5, and which have negligible effects-the accelerations in the system are everywhere much 



1 Since the term 'quasi- linear' partial differential equations is used by mathematicians in a different sense than 
here, I am now referring to such theories as 'practically linear' 
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larger than a , and hence we may put p = 1 everywhere within R. In such regions, the MOND 
equations coincide with the standard Poisson equation, and so the MOND potential satisfies the 
linear Poisson equation. However, it is not the standard, Newtonian solution, which assumes that 
the Poisson equation is valid everywhere to infinity. The fact that the field equation departs from 
Poisson beyond r M , imposes a different solution even at much smaller radii. 

It is easy to understand the effect with the aid of the so called 'phantom' mass (PM) distribution 
of the system, defined as (Milgrom 1986b) 

Pp= 4^G (A<A " A0Jv) ' (1) 
where <p and <p N are the MOND, and Newtonian, potentials of the system, respectively. In other 
words, the MOND correction to the Newtonian field is simply the Newtonian field of the PM. This 
concept is very useful in MOND generally, where it helps us bring to bear our long experience with 
Newtonian gravity, especially in cases where we have some prior knowledge of the PM distribution. 
The concept of PM was used effectively, e.g., in Milgrom (1986b), Milgrom & Sanders (2008), Wu 
& al. (2008), and in Zhao & Famaey (2010). It is especially useful in the context of QUMOND, 
where p v can be easily calculated at the outset. 

Quasi-staticity, in the sense defined above, has to hold for the auxiliary of PM to be useful, 
since we assume that its distribution is uniquely determined by the instantaneous distribution of 
p. Our nonrelativistic approximation breaks down if the system varies on time scales that are not 
much longer than r M /c. 

Under the above conditions, the MOND dynamics of any high-acceleration, spherical system 
are exactly Newtonian. This follows by applying the Gauss theorem to spherical volumes concentric 
with the system. The fact that asymptotically the field crosses to the MOND regime is immaterial 
inside the system. Viewed differently, for such a spherical system with R <C r M , the phantom mass 
forms a hollow, spherical cavity roughly beyond r M ; so it does not affect the dynamics of the masses 
within the system, which remain strictly Newtonian. 

When, however, the system p is aspherical, the distribution of p p is also aspherical. So, even 
if it is hollow inside r M , it can introduce finite MOND effects into the dynamics of our system. It 
is these effects that I want to calculate. 

Since R <C r M , the MOND field at, and beyond, r M is spherical to lowest order, with a 
correction of quadrupolar angular distribution. 2 Then the PM is also spherical to lowest order, 
plus a small aspherical contribution. The spherical part has no effect on the dynamics of the 
system. Since the aspherical PM lies beyond r M , its added potential much inside r M is of the 
quadrupole type (the dipole will be shown to vanish) 

K ~ CQ.jrr'. (2) 



2 I assume that the quadrupole of the system does not vanish. Otherwise, the dominant correction is of a higher 
multipole, and of a higher order in R/r M - 
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This is a vacuum solution of the Poisson equation since Qij is traceless. It is Qij that I want to 
calculate. 

To recapitulate, much inside r M the MOND acceleration field is very high; so there \i = 1 
(apart from the above mentioned, inconsequential, MOND islands). Thus, <p is, there, a solution of 
the Poisson equation. But it is not the standard solution, c/) N , which is oblivious of the MOND limit 
of the theory beyond r M . Instead, it is <f> ps <p N + <fi a , while at larger radii, r > r M , <\> becomes the 
full MOND solution of the problem. Because of the assumption of quasi-staticity, we can assume 
that even if the system changes over time, the distribution of PM 'adjusts itself continuously to 
the instantaneous configuration. 



2.1. Quasi-static approximation 

Even in deeply nonrelativistic systems, there are aspects that require relativistic treatment. 
Such is the case, e.g., when we consider influences over distances for which the light travel time is 
longer than time scales over which the influences vary. Such may be the case in the present context: 
The origin of the anomaly here is the behavior of the field at distances > r M , which, in turn, is 
determined by the matter distribution near the origin. If the latter changes on a time scale r that 
is not much longer than 2r M /c, the field at r M cannot be assumed to adjust to the instantaneous 
configuration, and instantaneously influence back the dynamics within p. A relativistic treatment 
is then needed, and my treatment below is not valid. 

In a system, such as a binary, whose p varies on a dynamical time scale r ~ R/v ~ (R? /MG) 1 ^ 2 , 
the condition for quasi-staticity can be written as 4(g N /a )(v/c) 2 <C 1. This can also be written as 

R » R R = {2R s r 2 M ) 1 ^ ~ (2vr J R2^) 1 /3 „ i(M/M ) 2 / 3 au (3) 

(R s is the Schwarzschild radius of M, and £h the Hubble radius). 

When the opposite of inequality (3) holds, one might be tempted to simply time average the 
distribution of the PM; but this has to be justified via a relativistic treatment. 



3. Calculation in QUMOND 

QUMOND (Milgrom 2010a) is a practically linear formulation of MOND derived from an 
action. It is the nonrelativistic limit of a certain formulation of bimetric MOND (BIMOND) 
(Milgrom 2009b). The field equation for the gravitational potential is 

= V ■ [v (^^j V<M, where = 4ttG P , (4) 

with V(f> — > at infinity. Here, v(y — > oo) — > 1, and v(y — > 0) — > y~ 1 ^ 2 . It is related to the 
usual MOND interpolating function, p(x), by u[xp(x)] = 1/ p(x). QUMOND thus requires solving 
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(twice) the linear Poisson equation (and not a nonlinear differential equation). It has already been 
put to good use for predicting and calculating MOND effects in the Solar system (Milgrom 2009a, 
Galianni & al. 2011), for calculating MOND fields of galaxies (e.g. Angus et al. 2012), and for 
structure formation in MOND (Llinares 2011); see also Zhao &: Famaey (2010). Its relation to the 
nonlinear- Poisson formulation is elaborated on in Milgrom (2012). For example, it is shown there 
that the two theories are equivalent for one-dimensional systems when the nonlinear Lagrangians of 
the two theories are the Legendre transforms of each other, leading to the above relation between 
p{x) and v{y). 

In regions where p = 0, the phantom density is given by 

and is easily calculated from the Newtonian potential. Since p p is nonzero only at distances much 
larger then R, we can take only the lowest contributing harmonic to the Newtonian potential. 
Taking the origin at the center of mass of the system, the dipole moment vanishes, and we keep 
only the monopole and quadrupole contributions, writing 

GM t 

4>n ~ 1- rj, (6) 

r 

where, 

V = (■'■ "rr'Q;,. Q i:j = IJ p(r')(r% - 3r#)dV '. (7) 

Using eqs.(6-7) in expression (5), keeping only terms up to first order in rj, we get: 3 

AttG Pp « -2yv'(y)^ + [18yi/(y) + 6yV'(y)]-J, (8) 

where, y = (r M /r) 2 . The first term is spherical. At large distances, where y <C 1, so u' — > — y~ 3 / 2 /2, 
it gives the asymptotic phantom matter 'isothermal halo': p p oc Mt/r M r 2 . At short distances, 
v 1 — > 0, and, as stated above, I assume that it does so fast enough that we can take it as zero within 
R. The first term thus describes a spherical distribution with a spherical cavity; so it has no effect 
on the dynamics of p. In the second term, call it 4irGp p , all the factors are spherical, except rj, 
which is a quadrupole. This part of p p also has an empty cavity for r r M , but it does produce 
a field inside it. Because p p is nonvanishing only at distances r > fi, we can keep only its lowest 
multipole contribution, which is the quadrupole, since rj is reflection symmetric. Its field within R 
can be written to the dominant power in r/r M as 

4> a « GrVQtf, where Q tj = ^J p P {v')(r l2 5 l3 - 3r#)^- (9) 



3 All the first-order terms are scalars that are linear in Qtj, so they must be of the form Gf{r)r l r^Qij = f{r)r 5 r] 
(Qu = 0). For example, use is made of r ■ Vrj — —3rj, and r % r 3, q,i,j = 12r?. 
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Inserting expression (8) for p p , with expression (7) for 77, we get 

9 f°° y f 

Qv = s^* Qkl J d vy 5/2 l u '(y) + 3 v "(y)] J dn ( 6 ij - 3n i n j) n k n h (10) 

where n = r'/r'. For the radial integral to converge at r — > 0, it is necessary that v\y) vanishes 
at large y faster than y~ 7 ^ 2 . This is, indeed, what is quantitatively meant when I say repeatedly 
that v is assumed to approach 1 fast enough. For slower vanishing, the MOND correction of type 1 
discussed in the introduction can be shown to dominate the present effect, and our approximation 
here is not valid. 

After performing the angular integral, and integrating the y integrals by parts, we get 

a 3 f°° 

Qij = — r Qij, where a = - dyy 3/2 [v(y) - 1]. (11) 

r M 4 JO 

Thus, a is a numerical factor that depends on the interpolating function. 4 

A generic choice of v(y) in MOND would be one for which v(y) — 1 is « y~ 1 ^ 2 for y <C 1, and 
can be made to vanish quickly for y S> 1. Such a choice would give a ~ 1, to within an order of 
magnitude roughly. 5 For example, for the very sharply transiting, limiting form of [i: \i{x) = x for 
x < 1, and fj,(x) = 1 for x > 1 [for which z/(y) = y^ 1 ^ 2 for y < 1, and v(y) = 1 for y > 1], we get 
a = 3/40. For z/(y) = (1 — e~ y )~ l l 2 one gets a ps 0.6. For the slowly transiting = (1 — e _s/1/2 )~ 1 , 
which was used successfully in Famaey & McGaugh (2012) for rotation-curve predictions, we have 
a 37. 

There may also be interest in interpolating functions for which very-near-Newtonian behavior 
is reached, not beyond g N ~ a , but beyond g N ~ na , with, possibly, k » 1 (see section 6). For 
these, v — 1 is made to vanish quickly only for while for y <C k, f — 1 ps y^ 1 ^ 2 - We can write 

some such functions in the form 

v = v K (y) = l + K- 1 ' 2 [v{y/K)-l], (12) 

where z/ is of the generic type defined above. Substituting expression (12) in eq.(ll), we see that 
for such a function a K = K 2 a K= i, where a K= \ is roughly of order 1. 

For a highly aspherical system, for which Q ~ MtR 2 (such as a system of two comparable 
masses), the anomalous acceleration is g a ~ GMtB? /r\ n and its strength relative to the Newtonian 
acceleration, g N ~ GM t /R 2 , is g a /gN ~ (R/r M ) 5 ~ (ao/yiv) 5 ^ 2 - This is so, I reemphasize, even if 



4 Note that Qij is the 'internal' quadrupole of the PM. It is finite, and picks up its main contribution from radii 
around r M - The expression for the 'external' quadrupole moment of the PM, defined as in eq.(7), diverges linearly 
at large radii (small y) for isolated systems, just as the total PM mass does. This does not lead to divergences in the 
fields, and is, anyhow, cut off by external fields. In some systems there are also departures from quasi-staticity to be 
reckoned with at large radii. 

5 Of course, a can be made arbitrarily large if v — 1 vanishes nearly as y -5 ^ 2 . 
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v {.9n/o>o) — 1 approaches zero arbitrarily fast with increasing g N /a . In almost spherical systems, 
the anomaly is further reduced. If m <C Mt is the 'aspherical' mass (e.g., the mass of planets in a 
planetary system) then the anomalous MOND acceleration in the system is g a ~ (m/M t )(R/r M ) s a . 

As an example, consider the MOND correction to the Newtonian, two-body force between 
masses m and M, a distance I apart, for which condition (3) for quasi-staticity holds. Here, 
Q zz = —£ 2 mM/M t , and Q xx = Q yy = —Q zz /2, where the masses lie on the z axis. The MOND 
correction to the Newtonian force is then attractive, and equals 

F a = 2a(£/r M ) 3 (mM) 2 M t - 3 a . (13) 



4. The nonlinear-Poisson theory 

In this theory, the MOND potential for an isolated system is the (unique) solution of the 
equation (Bekenstein Sz Milgrom 1984) 

V • [fi (^^J Vc/>] = 4nGp, (14) 

with V^> —7- at infinity. This is the nonrelativistic limit of Einstein-Aether theories (Zlosnik &: 
al. 2007), and it is also part of the nonrelativistic limit of TeVeS (Bekenstein 2004), and of generic 
formulations of BIMOND (Milgrom 2009b, 2010b). 

If p is related to u(y) of QUMOND as described above, the two theories coincide for spherical 
systems, and the (radial) dominant asymptotic behavior is the same in the two theories. Again, 
in the phantom-density approach, we need the dominant, asymptotic, aspherical MOND potential. 
However, in this theory we cannot calculate p p before a full solution of the problem is known. 

In Milgrom (1986a), I showed that the aspherical, far field still has, generically, a quadrupolar 
angular dependence, but its radial dependence is somewhat different from that in QUMOND. The 
asymptotic potential, beyond r M , is of the form 

0^(Ga o M t )VHn(r) + G^^, (15) 

where Sij is a symmetric, traceless (constant) matrix; it depends on the mass distribution p, but I 
do not know how. The phantom density outside p is 4irGp p = A<p, so its asymptotic form is 

iwGpp ^^ + fiG ^ (16) 

r M r A r (4+V3) 

where (3 = (v / 3 + 2)(v / 3 — 3). Here too, in high-acceleration regions, roughly within r M , p p vanishes. 
This is because there \i = 1, so Acp = A(j) N . The first term in the expression for the asymptotic p p 
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is spherical, and does not affect the dynamics in the system. The second term gives rise again to a 
quadrupole field 

' M 

S has the dimensions of mass x length^. It vanishes for a spherical system, and also in the 
limit a — > 0. So, we may write it as S ~ r M ), where R is the characteristic extent 

of the system. The dimensionless function £ (u) is unknown, however, and depends on the various 
dimensionless parameters of the system (mass ratios, ratios of distances, etc.), and has to vanish 
for u = 0, which corresponds to a = 0. Numerical calculations are needed to say more on this. 



5. Contribution of the small MOND domains near zero-gravity points 

Clearly, in every mass system there are critical points of the potential, where its gradient 
vanishes. 6 This follows from topological considerations, but is otherwise obvious. There are, for 
example, the points at (or very near) the centers of spherical stars or planets. Also, in a system 
of many compact objects, there are zero-gravity points, around which there is a MONDian region 
where \V<j>\ < a . These regions are surrounded by Newtonian regions, where we assume that 
p, = v = 1 to the desired accuracy. This implies that the total phantom mass within each such 
MONDian region vanishes: Applying the Gauss theorem to the volume V, within a surface X that 
is wholly in the Newtonian region, we have J v p p dV oc J E (V0 — V(p N ) • da = f s [p(\V <fr\ / a Q )V 4> — 
V4> N ] ■ da = J v V • 4>\ / a )V 4> — V<fi N ]dV = 0, where I used the fact that on S, p = 1. In a 
similar way it is seen that this is also true in QUMOND. 

The effects of such regions on the dynamics of p can be calculated if we know the PM dis- 
tribution in them. Since the typical acceleration within such a region is, by definition, |V</>| ~ a 
(I restrict myself to MOND theories that approach Newtonian dynamics quickly beyond a ), the 
characteristic size, L, and phantom density are related by AirGp p L ~ a . 

The MONDian regions near the centers of spherical objects are, themselves, spherical (they 
may be slightly aspherical if the general geometry is strongly aspherical, but I neglect this); so, 
since they have vanishing total PM, they have no outside effect. 7 

Other zero-gravity points typically occur in vacuum, so they cannot be extrema of the potential; 
they are thus saddle points. Take for example a system of a few comparable masses with inter-mass 
distances ~ R. If g(r) is the gravitational acceleration field in the system, g ~ GMt/R 2 , then at a 
zero gravity point, ro, we have g(ro) = 0. The size of the region around this point, within which 



I speak here of the nonrelativistic potential, therefore barring singularities, such as black holes. 

7 For a star of central density p, the characteristic phantom density in the region is also ~ p; therefore, the size of 
this region is L ~ a /Gp, which for p ~ lgr cm^ 3 , is a fraction of a centimetre. 
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g < a is L ~ a a (dg / dr)~ l ~ (a /g)R ~ R 3 /^. 8 The total PM in this region vanishes. To lowest 
order in R/r M the region is reflection symmetric; so its dipole moment vanishes to this order (it 
would, otherwise scale as p p L 4 ). The dipole moment is thus of the higher order: d ~ p p L 5 /R. The 
phantom quadrupole moment is Q ~ pL 5 . The acceleration field of the dipole at a characteristic 
distance i?, is ga ~ Gd/R 3 ~ Gp p L 5 /R 4 , and that of the quadrupole is g Q ~ GQ/R 4 ~ Gp p L 5 /R 4 . 
They thus scale in the same way. Since GLp p ~ a , we find that the anomalous MOND correction 
to the acceleration in the system is <7 a ~ a (L/R) 4 ~ a (R/r M ) 8 . This estimate can also be made 
for the case of a small mass m in a larger, Newtonian system whose characteristic size is R, and 
whose total mass is M t S> m. Again, one estimates that the MOND correction to the acceleration 
on m is g a ~ a (R/r M ) s . All this is confirmed by explicit calculations in QUMOND, which I omit 
here. 

The accelerations due to this effect thus scale as (R/r M ) 8 , compared with (R/r M ) 3 for the 
correction discussed above; but for this effect the anomalous acceleration on a small mass does not 
vanish in the limit m — > 0. 

6. Discussion 

Intuitively, one would expect that the MONDian, two-body force becomes exactly Newtonian 
at high accelerations, if the interpolating function does so. I have shown that this is not the case, 
and that, more generally, there are MOND effects that linger even in high-acceleration systems. 

For theories with a ~ 1, I cannot think, at present, of a precision experiment to test this effect: 
in the Solar system, the effect is too small to be tested with present capabilities. For example, for 
the dynamics of the massive planets, m/Mt ~ 10~ 3 , and r ~ R ~ 10 _3 r M ; so the added MOND 
accelerations are of order W~ 12 aa , where present accuracy hardly reaches ~ 10 -4 a . 

Note that my treatment here is not applicable, as is, to short-orbital-period pulsar binaries, 
as they are anything but quasi-static: For binaries that serve as laboratories to test relativistic 
theories, having orbital periods of a fraction of a day, the orbital light cylinder radius is of the 
order of a light hour, much smaller than the MOND radius of these systems, which is larger than 
10 3 light hours. For the Solar system, we see that condition (3) is satisfied by the large planets, 
which determine the quadrupole of the PM. 

But do we necessarily have a ~ 1? In considering MOND, I have always labored under the 
supposition that it involves only one acceleration constant that plays all possible roles of such 
a constant. This, to me, is a basic tenet of MOND. This supposition is amply supported by 
rotation-curve and other dynamical studies of galaxies, where indeed the same acceleration constant 



If some masses are much smaller than others, then the zero-gravity points are much nearer such masses than R; 
say a distance R m . Then L ~ (a /g)R m ; but our result below remains valid. 
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appears with different roles (see, e.g., Famaey & McGaugh 2012). 9 In particular, this means that 
a that appears in the mass-asymptotic-speed relation, MGa = V£, also marks the boundary of 
the Newtonian regime; namely, that beyond ~ a , Newtonian behavior is reached quickly. This 
underlies my conclusion above regarding the Solar system, since it implies that a ~ 1 within roughly 
an order of magnitude. 

This supposition holds (or can easily be made to hold) in most relativistic formulations of 
MOND: MOND adaptations of Einstein- Aether theories (Zlosnik & al. 2007), BIMOND (Milgrom 
2009b, 2010b), theories based on a polarizable medium (Blanchet & Le Tiec 2009), and nonlocal, 
metric theories (Deffayet k, al. 2011). However, it does not hold in TeVeS, in its original form. The 
basic reason behind this is that this form of TeVeS does not tend to general relativity (GR) in the 
formal limit a — > 0: it is not quite GR compatible. To get GR in such a TeVeS version, one also has 
to take a limit k — > 0, where A; is a dimensionless parameter of TeVeS (possibly representing several 
such parameters in different TeVeS formulations). Since various tests, in the Solar syste and pulsar 
binaries, strongly constrain high-acceleration departures from GR, they force one to take k <C 1. 
Now, in the nonrelativistic limit of such theories, the fully Newtonian regime can occur only beyond 
9n ~ a-l = K Oo) where k>1, Effectively, we then have two very different acceleration constants, a 
and al, playing important roles in the theory. For example, Bekenstein & Magueijo (2006) found 
that for the TeVeS formulation they used, the boundary of MOND regions around zero-gravity 
points are defined using k = (4ir/k) 2 S> 1 (k « 1.75 • 10 5 for their choice of k = 0.03). Starting 
from the same assumptions, Magueijo & Mozaffari 2012 used even higher values of k > 1.6 • 10 6 
(based on k < 0.01). Galianni &; al. 2011 follow a similar line in most of their paper, but do make 
the point that other, no less sensible, choices of a MOND theory, having k ~ 1, give effects that 
are many orders smaller. 

In my opinion, a strong desideratum of a relativistic MOND theory that involves only one 
dimensioned parameter, a , is GR compatibility: the theory should tend to GR in the formal limit 
a — > 0. This is analogous to the requirement that GR go to Newtonian dynamics for c — > oo; it is 
also analogous to Bohr's correspondence principle in quantum mechanics. All relativistic MOND 
theories mentioned above, except TeVeS, satisfy this requirement. Thus, in all these theories 
k ~ 1-which is preferred also by galaxy dynamics-is admissible. As already mentioned, the original 
versions of TeVeS are forced to adopt k> 1 because they are not GR compatible, a fact that puts 
them in danger of conflicts with tests in the high acceleration regime, unless k <C 1. Babichev & al. 
(2011) have recently devised a theory that might be viewed as a short-distance modification of TeVeS 
(introducing an additional scale length besides a ), and that avoids conflicts with observations, while 
admitting k ~ 1. I feel that a theory with k» 1, and all the consequences that follow from it, are 
not generic MOND results. 

The QUMOND versions of theories with k » 1 would involve an interpolating function as given 



This is similar to h playing all the roles of an action constant in quantum mechanics, or to c playing all the roles 
of a critical speed in relativity. 
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in eq.(12). Then, as we saw in section 3, the coefficient a defined in eq.(ll) is a = a K = n 2 a K= i, 
where a K= \ is of order 1, to within an order of magnitude, roughly. The high k values used in 
the above-mentioned analyses, could already be in conflict with existing Solar-system limits on the 
effect I discuss here. Checking this requires a detailed analysis. Here I give only a rough estimate: 
Assume that the quadrupole moment of the Solar system comes from a single, most dominant 
planet, of mass m and semi-major axis I. For convenience of comparison, write the anomalous 
potential as 

<p a = -^r l r> (eiej - j^). (18) 

Then, from all the above, we have 

Q 2 « -3K 2 a K=1 G£ 2 mr- 5 . (19) 

Taking the dominant contribution from Neptune (Jupiter, Saturn, and Uranus give only somewhat 
smaller contributions), we get 

Q 2 « -3 • W~ U K 2 a K=1 s- 2 . (20) 

In comparison, the external-field effect studied by Milgrom (2009a), and by Blanchet &; Novak 
(2011), produces Q 2 ~ (0.2 — 4) • 10~ 26 s~ 2 . So, if a K =i ~ 1, the effect discussed here is larger for 
k > 10 4 . Values of k > 10 5 might already be excluded. A more detailed analysis should consider 
that several planets contribute to the quadrupole of the Solar system, that these moments are 
rotating with the planets, that our results were derived in QUMOND, while the nonrelativistic 
limit of a theory like TeVeS is not QUMOND, that the Solar system is not isolated, but embedded 
in the Galactic field, as well as other factors. 
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